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ABSTRACT 

We examine the reflection and transmission of plane waves by a 
finite thickness shear layer having a linear velocity profile and 
bounded by two otherwise uniform parallel flows. The pressure per- 
turbation equation in the shear layer has been shown previously to 
have exact solutions in terms of Whittaker M-functions. It is shown 
that in addition to the angle of plane wave incidence and the relative 
Mach number of the flows bounding the shear layer, the scattering 
properties of the shear layer depend crucially on a parameter x in 
such a manner that the case x ^ 0 characterizes the long wavelength 
properties of the layer (with x = 0 being the vortex sheet) and the 
case X -»■ " characterizes the short wavelength properties of the layer 
(Geometrical Acoustics) . We have evaluated numerically the relevant 
Whittaker M-functions and using these values we have studied the 
behavior of the energy reflection coefficient for a number of 
cases (angle of incidence, Mach number and various values of the 
parameter, x) at which the corresponding vortex sheet is known to 
exhibit resonances (reflection coefficient -»■ °°) and/or Brewster 



angles (reflection coefficient 0) . We find that, unlike the vortex 
sheet, the finite thickness shear layer has no resonances or Brewster 
angles. Moreover we find that in general for x > 0.5 the amplified 
reflection regime degenerates into the total reflection regime of 
geometrical acoustics even in the cases in which the corresponding 
vortex sheet has resonances. In contrast for the region of ordinary 
reflection we find that in the cases in which the corresponding 
vortex sheet does not have a Brewster angle the values of the reflec- 
tion coefficient up to x = 1 follow quite closely those of the vortex 
sheet whereas for the cases for which the corresponding vortex sheet 
has a Brewster angle the magnitude of the reflection coefficient may 
be quite sensitive even to small changes of x in certain cases, 
contrary to some previous results of Graham and Graham^, The results 
of the present studies indicate that caution should be exercised in 
uncritically modeling a finite thickness shear layer by a correspond- 
ing vortex sheet, a practice usually followed currently in noise 
research and aeroacoustics. 



3 


a. . 

13 


b 

ff g 
i 
k 
k 

X 

iti 


P » P 


r 

sgn 

t 

w 

X, z 

A f H / C / D 
H 

M(z) 

M-functions 

Jl 


R = ± iR2 


T = + iT2 


R.P. 

U(z) 


NOMENCLATURE 
= (constant) sound speed 

= constant coefficients defined by eq. (12) 

= linear velocity profile slope [units; (time) 

= the two independent solutions (eq. (9)) of eq. (7) 

= wave vector of incident plane 
= x-component of wave vector ic 
3 

= second index of the Whittaker M-f unctions. 

= pressure perturbation 

= linear combinations of f and g. 

= position vector 
= "sign of" 

= time 

= z-component of velocity perturbation 
= coordinates 

= expressions defined by eq. (14) 

= heaviside function 
= local Mach number 

= upper fluid Mach number 

= Whittaker M-functions 
= Energy reflection coefficient 

= Complex reflection coefficient for the perturbation 
velocity potential. 

= complex transmission coefficient for the perturbation 
velocity potential, 

= "Real Part of" 

= local mean speed 
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= non-dimensional variable defined by eq. (5) 

= incident wave angle (see Fig. (1)) 

= sin0 

= non-dimensional parameter defined by eq. (8) 

= angular frequency of incident plane wave 

= perturbation velocity potentials for the lower dnd 
upper regions respectively. 
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I. Introduction 

The shear layer with a linear velocity profile has been the object 
of a number of studies. Kuchemann^ considered the stability of the 
boundary layer with a linear velocity profile. Pridmore-Brown^ 
studied rectangular duct modes in a duct with the basic flow having 
a linear velocity profile. Graham and Graham^ studied plane wave 
propagation through a linear veloicty profile shear layer. Goldstein 
and Rice** found an exact solution to the pressure perturbation equa- 
tion in terms of linear combinations of parabolic cylinder functions 
of different order. We have reported earlier (see Koutsoyannis®) 
some preliminary results on sound propagation through a shear layer 
with a linear velocity profile using essentially Whittaker M-functions 
as the basic solutions of the pressure perturbation equation. Jones® 
studied the stability of such a layer for subsonic basic flow using 
our previous solution in terms of Whittaker M-functions. Recently 
Scott^ examined wave propagation through a linear shear layer using 
the Goldstein and Rice** unwieldy solution. With the exception of 
Goldstein and Rice**, Koutsoyannis ® , Jones® and Scott^ the other 
investigators were concerned with either series or asymptotic solu- 
tions of the pressure perturbation equation in the shear layer region. 

In I the present study we are concerned with the behavior of the 
energy reflection coefficient for plane waves incident on a finite 
thickness shear alyer with a linear velocity profile and bounded by 
two uniform parallel flows. For that purpose we have evaluated 
numerically the solutions of the pressure perturbation equation in 
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terms of Whittaker M-functions and with the aid of these we have 
numerically evaluated the reflection coefficient for a number of 
typical cases involving the relevant parameters, i.e., the angle of 
incidence of the plane waves, the upper fluid Mach number and a 
characteristic parameter x which represents a non-dimensional measure 
of the disturbance Strouhal number with respect to the disturbance 
Mach number in the mean flow direction (see Section IV.)- 

In section II below we shall derive the expression for the 
reflection coefficient in terms of the two basic solutions of the pres- 
sure perturbation equation (i.e., the Whittaker M-function solutions). 
In section III we shall summarize and discuss our numerical results 
on the variation of the reflection coefficient. In section IV we 
shall discuss various limiting cases including the vortex sheet and 
geometrical acoustics limits and in section V we shall give an 
analytical proof for the absence of resonances for a thin shear 
layer. The results of the present studies are then summarized in 
the final section VI. At the end, in Appendix A we discuss the 
representations of the perturbations in the three flow regions 
(see Figure 1) , and the definition of the energy reflection coefficient 
as well as energy conservation and finally in Appendix B we give the 
general outline of the numerical scheme used to evaluate the Whittaker 
M-functions . 

II. The Expression for the Energy Reflection Coefficient 

Without loss of generality we may assume that the layer with the 
linear velocity profile is bounded by two uniform flows, one of which 
is at relative rest, as shown in Figure 1, i.e.. 


we assume the 
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following two-dimensional lin the (x,z) plane] inviscid compressible 
shear layer characterized by the mean continuous velocity flow field; 

U (z ) = 0 for z £ 0 

= bz for 0 £ z £ z, (1) 

= bz for z^ £ z 

Furthermore we assume a time-harmonic plane wave incident from the 

z ^ 0 half space with wave vector ic and wave number k = — . The entire 

a 

unperturbed flow field is otherwise assumed to be homogeneous, i.e. 
no variation for mean densities, temperatures and speeds of sound 
are allowed. 

The velocity potentials 4>o lower and upper regions 

iv XI 

of uniform flow, respectively, bounding the shear layer of thickness 

1 7 * 

z^ may be taken to be (see for instance Graham and Graham* or Miles'): 



exp 

±ik(xsin9 + zcos9 - at) 

+ 


( J 




(2) 

+R 

exp 

±ik(xsin9 + zcos9 - at) 

1 1 for z £ 0 


(p = R.P. < ±iT 

{exp 

|±ik(xsin9 ± (z-Zj^)|sin6 

'iV"! ■^■^'^1} 

1 for z > z 


* 

The representation of the perturbation fields expressed by equations 
(2), (3), (4) and (5) is essentially the one used by Graham and Graham^, 

although they did neither use the non-dimensional variable n nor the 
parameter t (see equations (5) and (8)). At any rate this representa- 
tion is consistent with the different ones used by Miles and Ribner® . 
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In the above equations R.P. denotes "real part of"; the first term in 
equation (2) represents the incident plane wave of unit amplitude 
emanating from the half space z ^ 0 having wave vector k that makes 
an angle 6 with the z-axis as shown in figure 1, with - ^ ^ 0 :< + -j; 

R and T are the complex reflection and transmission coefficients 
for the velocity potentials 4)^ and respectively, i.e.: 


R = ± i R 2 

T = ± i T 2 , 

where R^ and are real for the so-called ordinary and amplified 
reflection regimes and complex in the region of total reflection. 
In equations (2), (3) and (4) above the upper signs are 

to be used for ri^ > 1 (ordinary reflection regime) , and the lower 
signs for < -1 (amplified reflection regime); is the non- 

dimensional quantity ~ (see equation (5) below) with 


bz 


M, = — = 
1 a 


being the upper fluid Mach number. This choice of signs 


ensures that the radiation condition, as discussed by Miles 
Ribner® and Graham and Graham^ is satisfied. 


In the shear layer region 0 z ^ in which the mean flow 
varies linearly with z (see equation (1)), it is convenient to use 
the non-dimensional variable r\ (see Ref. 5) 


(4) 


n 


1 

sine 


bz 

a 


1 

sine 


M(z) 


( 5 ) 


where M(z) is the local Mach number of the mean flow in the shear 
layer. If then one assumes that within the shear layer region the 
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pressure perturbation p' (r,t) is of the form: 

^ (cos) 

p' (r,t) = p (z) sin (k^x-oJt) 

(0 

with k = k*e = — k sin0, one obtains the following ordinary differ- 
ential equation for the z-dependent part p(z) of the pressure pertur- 
bation p' in terms of the nondimensional variable n defined by 
equation (5) (see Koutsoyannis ^) : ' 

Pnn ■ f Pn * - Dp = o ,7, 

with t' being the parameter defined by: 

4t = ^ sine (8) 

We have shown earlier (see Koutsoyannis® ) that the two independent 
solutions of equation (7) are: 


f ( n ; T ) 
g ( n ; T ) 


(4iT)“ 


1 

n^M 


ix, ±m 


(4ixri^) 


m = I 


(9) 


which are real functions of the variable n and the parameter x for 
real values of n and x. Using the above solutions f(ri;x) g(ri;x) one 
may further write the following general expression for the pressure 
perturbation p' (r,t) and the z-component w of the velocity perturba- 
tion: 

p' (r,t) = p (n; x) sin jk (xsin0-at) I + p^^^ (ri; x) cosjk (xsin0-at) j 
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I P^^Nn;T)cos[k(xsine-at)] - (n; x.) sinjk (xsin0-at) ] | 

where p^^^(n) and p^^^(n) are linear combinations of the two inde- 
pendent solutions f (n) and g(n) (given by equations (9)) of the 
pressure perturbation equation (7), i.e. 

P^^Nn;x) = a^j^f(n;x) + , p^^Mn;x) = a2j^f(n;x) + a22g(n;x) 

with a. . constants. 

ID 

Finally we apply the appropriate boundary conditions at the two 
edges of the shear layer, at z = 0 and z = [or in terms of the 
non-dimensional variable ri (equation (5)) at ri^ = and 

Til = si^6 ' ” ^1^ ' namely continuity , of the pressure perturbation p' 

* 

and of the z-component w of the velocity perturbation. Then using 
equations (2), (3), (10), (11) and (12) one may obtain eight linear 

algebraic equations for the determination of the eight unknowns a^ j , 
and appearing in these equations. After a somewhat tedious 


Specifically at z = 0 we set: 

+ 

p' (z = 0 ) = -P^ = p' (z = 0 ) 

w(z = 0~) - (V4>^) . e^ = w(z = o'*’) 

and similarly at z = z^ we set: 

D(|> 

p'(z = Z^) = -p^ = p')Z . z^) 

^ 

w(z = z^) = (V(j)^) . e^ = w(z = z^) 

with D/Dt designating the convective name operator in the upper 
fluid region of constant Mach number M^. 


( 11 ) 


( 12 ) 
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but straightforward algebra we obtain the following expressions for 
the square of the reflection coefficient: 


R 


2 


(A±B)^ + (C±D)^ 
(A+B)2 + (C+D)2 


for I I ^ 1 



for I I £ 1 


(13) 


where the expressions for A, B, C and D given below involve essentially 
evaluation of the two independent solutions f(n) and g(n) of the 
pressure perturbation given by equations (9) and of their first 

derivatives at the two edges of the shear layer and are functions of 

11 * 

^o ~ sin¥ ' ’^1 “ sin^ ” ^1 parameter t (see equations (5), 

and (8) ) . 


* 

Using equations (4) for the definitions of the complex reflection and 
transmission coefficients R and T, respectively, we easily obtain: 


, for In^^l ^ 1 


i.e. in the ordinary and amplified reflection regimes where Rj^ and 
are real. In the total reflection regime | | ^ 1, Ri and Tj are 
complex and it is also easily obtained that in such as case |R|^ = 

I ^1 1 1 ^ 2 !^ T=0. The meaning of R^ follows from the inter- 

pretation: 


R = 


T = 


r2 + 
2 

+ 


"2 

,2 


_ Reflected Accustic Energy Flux (Time-averaged over one cycle) . 
~ Incident Acoustic Energy Flux (Time-averaged over one cycle) . 

and is discussed in Appendix A. 
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Explicitly : 


B = 4t ^Dq - 1 \f^ - 1 [f(l)g(0) - f(0)g(l)j 

C =yjr]^ - 1 [f(l)g^(0) - fn(0)g(l)j 

D = yJ l - T [f^(l)g(0) - f(0)g^(l)] 


(14 


The upper signs in equations (14) hold for > 1 and the lower 
signs for -If in both cases | I ^ If and we have used the 

notation 0 and 1 in the argvunents of f and g and their derivatives 
with the understanding that 0 designates evaluation at 
n = tIq = h I 2=0 “ 1 designates evaluation at n = = n | ^=2 ~ 

— i.e.f at the two edges of the shear layer. Equation (13) 


TT 


7T 


is valid for - ^ 6 ^ j with the upper signs holding for the regime 

2 

of ordinary reflection (g^^ > 1 and R < 1) and the lower signs for 

2 

the regime of the so-called amplified reflection (n^ ^ -If R ^ 1) . 
(For the total reflection regime |r| = If -1 ^ +lf 

- ^ 0 ^ + -j. ) It is seen from equations (13) and (14) that we 

recover the three reflection regimes; 


For simplicity we have assumed sin 0 ^ 0 in the expressions for A, B, 
C and D given by equations (14) since for the geometry chosen in 
Figure 1 resonances and/or Brewster angles for the corresponding 
vortex sheet cases exist only for 0^0^ Tf/2. Equations (14) may 
be made general f i.e. to apply for all values of the incidence angle 
0 (- ir/2 ^ 0 ^ + 7 t/ 2) by multiplying B and C by sgn(sin0). 
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2 

rij^ ^ 1 / R £ 1: Ordinary Reflection Regime 

2 

-1 £ £ t 1 , R =1: Total Reflection Regime 

2 

£ “1 / R £ Amplified Reflection Regime 

which agree with the results of Graham and Graham^. Indeed these 
regimes are the same as those found by Miles’ and Ribner® for the 
limiting case of the vortex sheet (t = 0) . 

2 

III. Numerical Evaluation of the Energy Reflection Coefficient R 

We have evaluated numerically the Whittaker M-functions involved 
in the two independent solutions f(n;T) and g(n;T) of the pressure 
perturbation equation given by equation (9) using the known series 
representations of f and g as given by Koutsoyannis ® for relatively 
small values of n and x whereas for large values of n and/or x we have 
used a numerical technique outlined in Appendix B. Using these values 
we have then evaluated the functions A, B, C and D given in equation 

2 

(14) and then the expression for the energy reflection coefficient R 
given by equation (13) . In particular we have chosen ranges of the 
relevant problem parameters, i.e., angle of incidence 6, upper fluid 
Mach nxmber M^^ and Strouhal number x, for which the corresponding 
vortex sheet (for the same 0 and M^) has one or two resonances and/or 
one or no Brewster angle. 

Figure 2 shows the behavior of the energy reflection coefficient 
R for the limiting cases of the vortex sheet (x = 0) and geometrical 
acoustics (x ■> «>) for a typical incidence angle 0 = 30°, as a function 
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of the variable r\-^ = = 2 - M^, or as a function of the 

upper fluid Mach number M^. As is known (see Refs. 7 and 8) the 

vortex sheet in this case has resonances at n, = - ^ and 

1 sine 

Tin ~ Q ,i.e. for M, = 4 and a Brewster angle at rii = + , 

1 COS0 1 ^1 COS0 

i.e. at =1.134. in the geometrical acoustics limit (t ^ ~) the 

2 

energy reflection coefficient R , as a function of ri^^/ degenerates ' 

• 2 2 
into the Heaviside-type step function R = 1 - H(1 -ti 2^)/ i.e. R =1 

2 

for < 1 and R = 0 for > 1. 

Figure 3 shows the variation of the reflection coefficient with 

upper fluid Mach number M^ for a niiitiber of values of the parameter t 

for a fixed angle of incidence of 30“ . For 30“ angle of incidence 

the corresponding vortex sheet (t = 0) has two resonances, i.e., 

at -1/sin 30“ and -1/cos 30“ and one Brewster angle at +l/cos 30“ . 

It is seen from the figure that even for very small t (-O.Ol) the 

resonances in the amplified reflection regime disappear and for 

values of t 0.5 the amplified reflection regime has degenerated 

2 

into the total reflection one with R = 1. In the ordinary reflection 
regime in which the corresponding vortex sheet has a Brewster angle 
at ri^ = +l/cos(30“) it is seen that even for low values of x ~ 0.05 
the Brewster angle has disappeared and there is discernible variation 
of the reflection coefficient with upper fluid Mach number up to the 
value T = 1 that we have calculated. 

Figure 4 shows the variation of the reflection coefficient with 
upper fluid Mach number M^^ for a number of values of the parameter x 
and for a fixed angle of incidence of 45“ . This is the special case 
in which the two resonances of the corresponding vortex sheet 
coalesce and the same holds for the Brewster angles. As in the 30“ 
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angle of incidence case we observe the same behavior in the amplified 
reflection regime > i.e., fast drop-off from the vortex sheet values 
with increasing x and by x > 0.5 the whole amplified reflection 
regime has coalesced with the total reflection one. In the 
ordinary reflection regime though we find essentially no variation 
from the vortex sheet values up to the value x = 1 calculated. 
Apparently it takes quite substantial values of x > 1 for the reflec- 
tion coefficient to start differing from the values of the correspond- 
ing vortex sheet. Finally Figure 5 shows the variation of the 
reflection coefficient with upper fluid Mach number with x as the 
varying parameter and for a fixed incidence angle 9 for which 
sin9 = 0.8. The corresponding vortex sheet has only one resonance 
and no Brewster angle. The general behavior is similar to that 
observed in Figure 4, i.e., rapid decrease of reflection coefficient 
values with increasing x in the amplified reflection regime and 
imperceptible change from the vortex sheet values in the ordinary 
reflection regime up to the value x = 1 calculated. 

From the above figures the following general conclusions may 
be drawn. The finite thickness shear layer shows no resonances and 
no Brewster angles. Moreover, the values of the reflection coeffi- 
cient in the amplified reflection regime drop-off substantially 
from those of the corresponding vortex sheet with increasing x, even 
for modest values of x on the order of ~1, whereas in the ordinary 
reflection regimes up to x ~ 1 modest or no significant changes are 
observed in the reflection coefficient values for the cases in which 
the corresponding vortex sheet exhibits no Brewster angle. In con- 
trast when the corresponding vortex sheet exhibits a Brewster angle 
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in the ordinary reflection regime the reflection (and transmission) 
characteristics of the finite thickness shear layer, as may be 
surmised from Figure 3, strongly depend on the value of the para- 
meter T (in addition to the angle of incidence 0 and the upper fluid 
Mach number M^) . This is in contrast to the results of Graham and 
Graham^ who for the case of = 3 and sin0 = 0.2, that they examined 
as a typical one for the ordinary reflection regime, they have calcu- 
lated negligible effects due to the finite layer thickness as compared 
with the corresponding vortex sheet results. The reason is, as may 
be easily checked, that in their case the corresponding vortex sheet 
has Brewster angle at n = =1.1021, i.e. very close to 1 (see Fig. 

2) and moreover they have evaluated the variation of R and T as a 
function of the non-dimensional parameter*; 

w , 

Shear Layer Thickness _ a ^ ^ ^ cos0 

Wave Length 2n b ‘ a ' 2n 

_ 0) „ COS0 
■ b ^1 2 tt 

and they considered variations of this parameter from zero to less 
than 0.2. This range corresponds to variation of our non-dimensional 
parameter x from zero to less than 0.02 and since the value x = 0 
characterizes the vortex sheet, it follows that the choice on the 


* 

In actuality the parameter in Ref. 1 is: 

Shear Layer Thickness _ 1 

(Incident wave vector component in the z-direction) ^ w cos0)”^ 
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above parameter in Ref. 1 has been an unfortunate one, and the 
range of values of this parameter investigated by Graham and Graham^ 
is too close to the vortex sheet value x = 0. 


IV. Limiting Behavior of the Reflection Coefficient 


In this section we study analytically the limiting values and 

2 

forms of the energy reflection coefficient R , as given by equation 
(13), as the parameters x, 6 and take extreme values by examining 
the behavior of A, B, C and D in equations (14) for these limiting 
values. In particular we will be interested in x ->■ 0 or x -»■ «> 
corresponding to the vortex sheet and geometrical acoustics limits, 
respectively (see Koutsoyannis ^) , in 0 0 or 0 -► ± tt/2 corresponding 

to the normal and parallel incidence, respectively, and in ->■ 0 or 
-»■ 00 , i.e., the low and high upper fluid Mach number limits, 
respectively. 

The vortex sheet limit x 0 follows from the known properties 
of the series solutions f (n) and g(ri) in equation (9) and from 
equations (14) (see Ref. 5). In this limit A 0 and B -»■ 0, whereas 


3n 


2 K sgn K 


/I - 



•k 

Since the parameter used in Ref.J^ does not contain the shear layer 
profile slope b, one could have anticipated that it is not a suitable 
non-dimensional parameter that adequately characterizes the scattering 
characteristics of the shear layer. In contrast our non-dimensional 
parameter x as given by equation (8) does characterize the essential 
features of the problem since it is the product of the Strouhal num- 
ber w/b and the sine of the angle of the incident wave vector with 
the z-axis. (See also Ref. 5 for a detailed discussion on the' physical 
meaning and importance of the parameter x.) 
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and 

D ^ -3n^ , 

and since k = sinS, obtain from equation (13) 


C ± D 


C + D 



(1 - Mj^sinS)^ ± sec0 ^(1 - M^sinS)^ - sin^0 | (15) 


(1 - Mj^sinB) + sec0 


^(1 - M^si 


in0 ) ^ - sin^0 


where again is the upper fluid Mach number and upper signs corre- 
spond to ordinary and lower signs to amplified reflection. Equation 
(15) is the equation for the vortex sheet and agrees with the corre- 
sponding results obtained by Miles and Ribner® . 

The limiting form of the reflection coefficient in the limit of 
(geometrical acoustics) may be obtained using the asymptotic 
form of the solutions f(ri) and g(n) in equation (9) as t ->■ “ (see 
Ref. 5). Instead of using these forms we may argue as follows: 

Since it has been shown that in the limit t ->■ «> one recovers both 

the amplitude and the phase function of geometrical acoustics (see 

2 

Ref. 5) it follows that the reflection coefficient R in the limit 
T 00 has the form consistent with geometrical acoustics, i.e., 


R^ = 0 for > 1 

= 1 for n-j^ < 1 


(16) 


The limiting form for normal incidence, i.e., 0 = 0 is obtained 
by observing that in this limit 


*^0 sin0 


"^1 sin0 


- Mj^ -> + 00 


-> + 00 


r 
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and also in equations (14) 

A -»■ 0 , B 0 and C ->■ -D. 

Consequently since only ordinary reflection is possible > 1) the 
reflection coefficient is zero. 

The limiting form for parallel incidence 6 = ± tt/2 i.e., 

2 

sin 6 = 1 is obtained from equations (14) by observing that in 
this limit 


A -> 00 and C ->• oo 
and the reflection coefficient 


^22 
2 A + C 

^ 2 2 
A + C 


The limiting form for M ^ 0 is obtained again from equations 
(14) by observing that in this limit 


A ^ 0, B -»■ 0 and C ^ -D. 

Consequently again only ordinary reflection is possible 

(n, n = — > 1) with the reflection coefficient being zero. 

' '1 'o sinS ^ 

Finally the limiting form of the reflection coefficient for 

M ^ 00 is the same as that for the vortex sheet, equation (15) , i.e., 
2 

R = 1 and since as M oo - oo only the amplified reflection 

regime applies. 

We summarize the above limiting forms of the values of the 
reflection coefficient together with the corresponding regimes 
(Ordinary, Total or Amplified reflection) in the following table: 
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T 

9 

Ml 

0 

00 

0 

± tt/2 

a 

00 

Reflection ^ 

Coefficient R'^ 

Vortex 
Sheet Limit 
Eq. (14) 

Geometrical 

Acoustics 

Limit 

0 

1 

0 

1 

Reflection Regime 

0. , T. , A. 

0. , T. , A. 

0. 

T. A. 

0. 

A. 


In the above table in the bottom line 0., T. or A. stand, respectively, 
for Ordinary, Total or Amplified reflection. 

V. Resonances for the Finite Thickness Thin Shear Layer 

In this section we present a proof for the nonexistence of 

resonances for a nonzero but thin shear layer. It is seen from 

2 

equation (13) for the reflection coefficient R that resonances 
which may only exist in the amplifying reflection regime < ~1) 
imply that 

A + B = 0 and C + D = 0 (17) 

We found out in deriving equation (15) for the vortex sheet case 
(t = 0) that, to the lowest order inx,A->-0 B ->0 whereas C and D 
are nonzero and yield the vortex sheet result. We next evaluate A 
and B to lowest order in t as follows; We first insert in equations 
(14) the values for the functions f(ri) and g(ri) from equation (9) to 
the lowest order in t and then we evaluate A + B at the zeros of 
C + D and we show that the two equations (17) are incompatible for 
small but finite t, the only possibility remaining being that of 
T = 0, i.e. the vortex sheet case. 
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For small t one may obtain from equation (9) (see Ref. 5) : 

f(ri;T) = 1 - + o(T^) 

o j/5 7\ , (18) 

g(n;T) = (to ' h) ^ ) 

Inserting these values in equation (14) we obtain: 

A = (4t) [3(riQ - n^)nQnj^(l - nQn^)j '' 

B. (4X,V^ (n^ nj) 

Observing that is the upper fluid Mach number we see 

from equation (19) that the condition A + B = 0 yields: 

A + B = (4X) |3riQnj^(l - ngHj^) +\/r\l - 1 - 1 (hq + nj + non;^) ] = ° 

Equation (20) must be evaluated at the vortex sheet values of 
i.e., where C + D=0, i.e. at the vortex sheet values of and ri^^* 
Using the two known properties of the vortex sheet solution, i.e. 

Oq + Hi = , nQOjL = 1 - ^M^ + 1 

Equation (20) yields: 

A + B = - (4t) (2M^) = 0 (21) 

We see that equation (21) is satisfied only if t = 0 or = 0 or 
both and consequently the two equations (17) required for the 
existence of resonances are incompatible to the lowest order in t ; 
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it follows that the nonzero thickness thin shear layer has no resonances 

* 

and no Brewster angles. 

One may also prove, in general, on the basis of certain differential 
properties of the quantities A, B, C and D, equaitons (14) , and the 
Wronskian of the solutions f(n;T) and g(n;T) (equation (9)) of equation 
(8) , that not only the thin shear layer, but also the finite thickness 
shear layer has no resonances or Brewster angles andmore over resonances 
and Brewster angles are possible if and only if either x = 0 (the 
vortex sheet case) or x -»■ “ (the geometrical acoustics limit (see. 
figure 2 ) .* ** 


VI. Conclusions 

We have evaluated numerically the energy reflection coefficient 
for plane waves incident on a plane shear layer having a linear 
velocity profile. The numerical computations were based on a 


* 2 
The conditxons for the existence of Brewster angles, i.e. R = 0, 

is the same as that for the existence for resonances, i.e. equations 

(17) , except that in this case the ordinary reflection regime (rii > 1) 

is the relevant one and the proof follows as for the case of the 

resonances. 

** 

This is true not only for the linear shear layer (equation 1) but 
also for a finite thickness shear layer of a generally continuous 
velocity profile. This result and the corresponding studies will 
be published in a separate following publication concerned with the 
scattering characteristics of finite thickness shear layers with a 
continuous velocity profile a special case of which is the shear 
layer with a linear mean velocity profile of the present study. 
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representation of the pressure perturbations in the shear layer region 
in terms of Whittaker M- functions. 

We have found that the shear layer exhibits no resonances and no 
Brewster angles and a separate analytical proof for the nonzero 
thickness but thin shear layer substantiates the absence of both 
resonances and Brewster angles for a finite thickness layer. More- 
over we have observed that, the behavior of the reflection coeffi- 
cient depends crucially on the parameter t, which ,as is seen from 
equation (8) , represents a non-dimensional measure of the disturbance 
Strouhal number with respect to the disturbance Mach number in the 
mean flow direction. In particular for moderate values of t the 
amplified reflection regime degenerates into the total reflection 
one whereas in the ordinary reflection regime the variation of the 
reflection coefficient with x depends on whether the corresponding 
vortex sheet has a Brewster angle or not. In cases in which the 
corresponding vortex sheet has a Brewster angle the reflection 
coefficient is sensitive to changes in x even for moderate values of 
X whereas in the cases in which the corresponding vortex sheet has 
no Brewster angle the reflection coefficient for moderate values of 
T follows rather closely the corresponding vortex sheet values. 

The above results indicate that caution should be exercised in 
modeling planar shear layers by vortex sheets uncritically even in 
the ordinary reflection regime and even for subsonic relative flows 
of the two regions bounding the shear layer, a practice customarily 
followed in current research and applications in noise studies and 
in aeroacoustics in general (see Ref. 9) . Although it is well known 
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that /the directional characteristics of the reflected and transmitted 
plane waves scattered by a finite thickness plane parallel 
shear layer are independent of the details of the velocity profile 
in the layer (provided that the profile is "smooth") (see for instance 
Ref. 10) , the amplitudes of the transmitted and reflected waves 
crucially depend on the parameter t, in addition to the angle of 
incidence 0 and the relative Mach niamber M of the two uniform flows 
bounding the finite thickness shear layer. Indeed we have shown that 
for certain combinations of t, 6 and M, i.e. those for which the 
corresponding vortex sheet has a Brewster angle (energy reflection 
coefficient = 0) the scattering characteristics of a finite thickness 
shear layer may drastically differ from those of the corresponding 
vortex sheet (t = 0) even for modest (non-zero) values of t. Moreover 
if one allows for different densities and temperatures in the three 
flow regions (z i 0, 0 ^ z ^ z^ and z >. z^^), a case of particular 
interest to noise generation by and/or propogation through hot jets, 
the difference between the scattering characteristics of a finite 
thickness shear layer and the corresponding vortex sheet characteristics 
becomes even more pronounced for certain ranges of values and/or 
combination of the relevant parameters involved. The detailed calcu- 
lations and the corresponding studies will be published in a following 
separate publication concerned with the scattering characteristics of 
a linear shear layer as examined in the present studies but allowing 
for different densities p and speeds of sound in the three regions 
of mean flow (see Figure 1) . 
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APPENDIX A 


A. 1 Representation of the Velocity Potentials in the Regions 

z £ 0 and z ^ z^. 

(i) Velocity Potential for the Lower Region z £ 0. 

Consider an incident plane wave which in coordinates 
(x,z) fixed in the lower fluid, which is at rest, is a sine 
wave with wavevector k making an angle 6 with the z-axis, i.e. 

ic.e = kcos0 
z 

with 1 9 1 • 

The wave is an upcoming one with amplitude A and thus its 
potential may be written as: 



incident 


A sin (k X +k z - wt) . 
X z 


And since k = ksin0 , k = 

X y 

we may further write: 

(|>^ = A sin |^(xsin0 


kcos0 , k = — 

Si 


+ ZCOS0 - at) 


The reflected wave will in general have in phase and out of 
phase case components; thus since the x-wave number is conserved 
in this stratified medium we may write for the reflected wave 
which is a downgoing wave: 


*^reflected ^ 


R^sin — (xsin0 - zcos0 - at) 

X Si 


+ R^cos — (xsin0 - ZCOS0 - at 
2 a 


A- 2 


where and R 2 are respectively the in-phase and out-of-phase 
components of the reflection coefficient for the velocity 
potential. Finally 

(J)p = Total Velocity Potential in the lower region z £ 0 


({).+(}) 

r 


= A I sin |^(xsin6 + zcosQ - at) j + 
+ Rj^sin |^(xsin6 - zcos6 - at) | + 
+ R 2 C 0 S |^(xsin0 - ZCOS0 - at) | 


Defining by R = R^ ± i R 2 » the complex reflection coefficient 
and taking the amplitude A to be real we may write in complex 
notation; 


1“ 

±i (k X 

e X 

+ k z 
z 

- wt) 

+ 

5 ± (k X - 
Re X 

k z - 
z 

.t,|| 

1“ 

GOS(k X 

+ k^z 

- wt) 

+ 

isin (k X 

+ k z 

- wt)] 

1 

X 

z 


X 

z 

iRo) 

cos (k X 

- k z 

- wt) 

+ 

isin (k X 

- ■ k z 

- <0t)]| 

z 

X 

z 


X 

z 


A I si 


sin(k X + k z - wt) + R, sin (k x - k z - wt) + 
X z l x z 


I 


+ R-sin (k X - k z - wt), 
2 X z j 


irrespective of upper or lower signs. 

(ii) Velocity Potential in the Upper Region z £ z^. 
The upper medium is moving with uniform velocity 
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U, = bz, ex, with = — — , with respect to the (x,z) - 

J. X J. a 

coordinate system which is fixed in the lower fluid (which 
is at rest) . 

Since again the x-wave number is conserved in the 
stratified medium we may write for the velocity potential 
of the transmitted wave: 

‘*’u jTj^sin (z-z^) - u)t)j + 

[^Tx^ 

r 

where and T 2 are the in-phase and out-of-phase components 

/ 

of the transmitted wave. As before is obtained from x-wave 

number conservation, i.e.: 

^Tx ^ ^ix ^ ksine = ^ sine, 

where ^ is the transmitted wave angle, but the z-wave number 
k^^ of the transmitted wave may only be uniquely specified by 
applying the radiation condition as postulated by Miles [17]. 
Namely from geometrical considerations. 



k^ is then obtained from conservation of the frequency w, 
i.e.f since^ is the angle that the transmitted wave vector 
makes with the z-axis, then 

0 ) = ka = k^c^ = k,j,(a + U^sin^) , 
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which together with the x-wavenumber conservation, and 
phase speed conservation, 

-^-rr = ^ + U, , 

sin0 sin 1 

results in: 


= k (l-Mj^sinS ) . 


Finally: 

(1-M^sin6)^ - k^sin^S 
= ± k I sin6 I - 1 . 

Which of the two signs in the expression for k^^ 

to be taken, is determined by the radiation conditions as 
postulated by Miles [7] i.e. + sign for > 1, i.e. in the 
so-called ordinary reflection regime and - sign for < -1 
i.e. in the so-called amplifying reflection regime. 

Finally writing 



T 



i T 


2 


and substituting in the expression for we obtain the 
expression given by eq. ( 3 ) in terms of the complex trans- 
misstion coefficient T. 
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A. 2. Definition of the Energy Reflection Coefficient. 

2 

We define the energy reflection coefficient R as 


2 _ Reflected Acoustic Energy Flux Density Averaged Over a Cycle 
~ Incident Acoustic Energy Flux Density Averaged over a Cycle 


i 

And since for a plane wave propagating in a medium at rest 
the energy flux density is aPQv' , where v is the (acoustic) particle 
velocity in the direction of the plane wave propagation, we may write: 


v-2 

incident 




2 


V 


,2 

reflected = 




thus indicating with a bar time averaging over a cycle of the incident 
monochromatic wave we may further write: 


V ' 2 

_ refl. 
incid. 



cos^ ^^(xsin9 + ZCOS0 - at)j + R^ sin^ | ^(xsin6 - zcos0 - at) j ^ 




+ R 


2 

2 * 
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A. 3 Acoustic Energy Conservation . 

Much confusion has resulted in the literature concerning the 
precise expression and the meaning of acoustic energy conservation in 
a stratified fluid since the original publication of Blokhintzev' s 
work [11]. The confusion has been compounded by Ribner's [8] 
analysis as well as by the more recent review by Candel [12] . It 
is actually a very simple matter to show that in parallel flows what 
is actually conserved is the time-averaged (Over one cycle) cross-flow 

component of the energy flux density in a reference frame moving with 

0 

the mean local speed, i.e. the conservation principle may be stated 
as follows: 


p ' v' . n^_ 

— ^ = const, (independant of z)', 

where (o' accounts for the usual Doppler factor relating (O (in the 
(x-z)-frame) and (o' (in the moving frame) and iij_ is the unit 
normal to the parallel flow direction. The above relation relates 
and compliments Bretherton and Garret 's[ 13] action principle with the 
essential difference that whereas the above principle is limited to 
short wavelengths, our result in the above equation. although limited to a 
statement concerning the time-averaged cross-flow part of the energy 
density flux of a single monochromatic component, it is valid for 
all frequencies. 

Since the proof of the above statement has hot been given any- 
where in the literature, we present below a sketch of this proof, 
with the details to be published elsewhere. 
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In a medium at rest 


P’ 



and v' 


y<)> , and 


p ' v' 



d(j> d({) 

^0 3t ^ ' 


thus for the lower region (which is at rest) we may write, since w = w' , 


p' v' . e^ 


W ' /(JJ 


-Pq 34.3;, 34 >J^ 


0) 


■5T~ T: 


2 

2^ COS0 (1-R ) . 


Similarly for the upper region of uniform flow = hz^, for which 


Di<f) 




l'"u 

XT • •= 


0 Dt 


, v' = V4> 


u 


we obtain similarly: 


~Pq ^ 

a)'/w / u' Dt 3z 

4i 


A^Pq wd-Mj^sinS) 


0 ) ' 


'‘tz 


2 



Thus energy conservation follows. 


a-nh k.^ = k^^ 


IZ 


where 


kcos0 
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APPENDIX B 


Numberical Evaluation of the Solutions of y" + 




y = 0 


We are concerned with the numberical evaluation of the solutions 
of the differential equation 

2 


y" + 


1 T . 1 - 
4 ■ X 

4x 


y = 0 


(B-1) 


subject to the coniditon 


1+y 

1 im y { X ) /x =1 

x-»-o 


(B-2) 


for arbitrary values of x > 0, i, and y 7 ^ -1, -2, -3,... 

Although power series and asymptotic series exist, they are 
not really useful for numberical work, especially for "large" values 
of T and x. A much better approach is to use the three term recurrence 
relation and normalization sum for numerical work. This approach 
allows a numerically stable, unitary, and computationally efficient * 
algorithm which can be used over a very wide range of values of x and T. 

The solutions of equation (B-1) and under dondition (B-2) are 

ity 

S(T, + y;x) = ^ 

1-y 

S(T,.- y;x) = M^x,-y/2 ^ 

where M, ^ (z) is Whittaker's function. 
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Using Buchholz [14] we can show that the functions 

u. = s(T,v + 2jfx) }:j;; i 
satisfy the three term recurrence relation 
y + j - 1 


(y + 2j - 1) (y + 2j - 2) 


1 + 
X 




2t 


(y + 2 j - 1) (y + 2 j - 2) 


u 


■ 

j(y + 2j 

2 2 
+ 1) + T 

(y + : 

j) (y + 2 j 

+ 1) (y + 2j + 2) 


= 0 


J^j-1 


(B-4) 


and are subject to the normalization relation 


(w) u . 
j = 0 J J 


wx 

2 


Ity 

2 
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The Aj (u) are given by 




Aj (w) = (i) J Pj 2 


(-iw) 


(B-6) 


The formulas (B-4) , (B-5) , (B-6) form the basis of the numerical 

algorithm. Using the Perron-Kreuser theorem [15] , the solutions of 
(B-4) behave asymptotically as 




Using (B-5) and (B-6 ) , we infer that the required solution has the 
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latter behavior; this solution is called the subdominant solution. 
The actual numerical algorithm is based upon the algorithm of Olver 
[16], 
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